Dendritic to globular morphology transition in ternary alloy solidification 
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The evolution of solidification microstructures in ternary metallic alloys is investigated by adap- 
tive finite element simulations of a general multicomponent phase-field model. A morphological 
transition from dendritic to globular growth is found by varying the alloy composition at a fixed 
undercooling. The dependence of the growth velocity and of the impurity segregation in the solid 
phase on the composition is analyzed and indicates a smooth type of transition between the dendritic 
and globular growth structures. 
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Multicomponent alloys form the most important class 
of metallic materials for technical and industrial pro- 
cesses. Combined with the number of components is 
a wealth of different phases, solidification processes and 
pattern formations. The mechanical properties of a mate- 
rial strongly depend on the morphology and on the char- 
acteristical quantities of the microstructure. One of the 
most common growth morphologies in metallic alloys are 
dendrites. During growth, dendritic patterns generate 
solute microsegregation forming the structure of grain 
boundaries on a larger scale. The growth of dendrites 
from an initially undercooled metallic melt is caused by 
a negative concentration gradient ahead of the solidifica- 
tion front. Extensive studies have been made to explore 
dendritic growth in pure substances and in binary alloys. 
The constitutional undercooling and the corresponding 
concentration gradient at the growing front are the main 
factors leading to morphological instabilities of the solid- 
liquid interface 0, 01 • In a binary alloy, the relation be- 
tween the temperature interval ATq of the liquidus and 
solidus line, the concentration of the impurity compo- 
nent Co in the melt, the partition coefficient fee and the 
liquidus slope is given by ATq — — meCo(l — ke)/ke, 
according to the equilibrium phase diagram of a station- 
ary planar front. The temperature interval determines 
the magnitude of the constitutional undercooling at a 
growing solid-liquid interface and is hence a significant 
material quantity for the stability properties. The more 
pronounced concentration boundary layer of the impurity 
ahead of the moving interface at larger values ATq drives 
the occurrence of instabilities leading to the formation of 
cellular and/or dendritic structures. For multicomponent 
alloys as e.g. in ternary systems, the interface instabili- 
ties are driven by the characteristic diffusion properties 
of each solute field in the multicomponent system. Us- 
ing an approximated model with independently acting 
solutes, it was analytically shown in that each so- 
lute enhances morphological instability of a solidification 
front with different weighting factors. A variation of the 
alloy composition causes a change of the strength of insta- 
bility forces on the growing interface, because the value 
of the liquidus-solidus temperature interval strongly de- 



pends on the alloy composition. Investigations of the in- 
teracting quantities and of the fundamental mechanisms 
during solidification of ternary or higher component ma- 
terials have just recently become more intensive by means 
of numerical modeling. 

To numerically treat the complex solid-liquid interface 
geometry of dendritic crystals, a phase-field approach has 
the advantage of avoiding the explicit tracking of the 
phase boundary. The solid and liquid are distinguished 
by a phase-field variable (j){x, t) with a smooth transition 
from to 1 leading to a diffuse interface. In the sharp 
interface limit, classical free boundary problems and the 
Gibbs-Thomson law can be recovered from the phase- 
field approach (see e.g. Q and references therein). In 
the two past decades, there has been a great progress in 
phase-field modeling of growth structures in pure sub- 
stances and alloy systems. For a historical background 
and applications see the review article in 5] and refer- 
ences therein. 

In this Letter, the influence of changing the alloy com- 
position on the interface stability, on the characteristic 
morphology and on the growth velocity is investigated by 
performing numerical simulations of a general multicom- 
ponent phase-field model with interacting and coupled 
diffusion fields. We chose the ternary Ni6oCu4o-a;Cr2, al- 
loy system as a prototype for this study and herewith 
build upon intensive studies of the binary Ni-Cu sys- 
tem (e.g. and 0)- Hence, the corresponding phys- 
ical parameters are established relatively well. Further- 
more, the liquidus-solidus interval for binary Ni6oCu4o 
is about three times larger than the interval for binary 
NigoCr4o and one can expect a significant dependence 
of growth morphology on the concentration of the alloy 
components. 

For the subsequent computations of ternary alloy solid- 
ification, we apply a recently developed phase-field model 
to the case of 3 components and 2 phases. The general 
formulation of the phase-field model in allows for an 
arbitrary number of phases and components and is based 
on an entropy density functional of the form 
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The bulk entropy density s depends on the internal en- 
ergy e, on a concentration vector c — (ci, . . . ,cjv) and 
on a phase-field vector <j) — The vector 

components c^, i = 1,...,N, represent the solute con- 
centrations in an iV-component alloy. The variable 0ct, 
a — 1, ... , M, denotes the local volume fraction of phase 
a in an M-phase system. The thermodynamics of the 
interfaces is determined by the gradient energy density 
a{(j),\7(j)), by the multiwell potential w{(f>) and by the 
small scale parameter e related to the thickness of in- 
terface. The gradient energy and the multiwell potential 
depend on the surface entropy density 7. The evolution 
equations are obtained from Eq. in a thermodynam- 
ically consistent way ensuring a combination of positive 
local entropy production and mass conservation laws. By 

variational derivatives dt(t)a ^ A, the phase- field 

equations read 

/ da da \ 1 dw I df 
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where /(c, T, (p) is the free energy density, A is a Lagrange 
multiplier taking the constraint (f)i + . . . + — 1 for the 
phase fields into account and a; is a kinetic factor related 
to the kinetic coefficient ^. In the sharp interface limit, 
it has been discussed in |^ that Eq. (0) relates to the 
Gibbs-Thomson equation for a moving interface. 

In our numerical investigations, the alloy composi- 
tion will be changed by keeping the Ni concentration 
at a constant 60at.% and by exchanging Cu by Cr, 
i.e. Ni6oCu4o-a;Cra; with < x < 40at.%. Three as- 
sumptions are applied: (i) We examine primary dendritic 
growth involving only two phases in the ternary system: a 
solid and a liquid. The eutectic region in the Ni-Cr phase 
diagram is not considered. For a solid-liquid system, the 
phase-field model reduces to one variable cj) denoting the 
fraction of solid phase, (ii) The system is considered in 
ideal solution approximation with a free energy density 
of the form 
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where Li and Ti are the latent heat and the melting tem- 
perature of component i, respectively, R is the gas con- 
stant, Vm is the molar volume and h{(j)) is a monotone 
function on [0, Ij that satisfies /i(0) = and h{\) = 1. 
(iii) Following |g , we impose an isothermal temperature 
field T . This approximation is valid because the thermal 
diffusivity of the alloy is about four orders of magnitude 
larger than the solutal diffusivity. In addition, for the 
considered undercooling, the growth takes place under 
diffusion limited conditions. 

Surface energy and kinetic anisotropy are incorporated 
into the model by formulating orientation dependent gra- 
dient energy a((/), V0) and kinetic factor V0) 

a(0,V0) =7|V0|2 = 7o(l-|-eeCOs(40))|V(/.p (4) 



t^(0,V(/)) = t^o (1 - Cfc cos(46l)) , (5) 

where 6 is the angle between the vector V(/) and the x 
axis. 

Assuming the mass fluxes to be linear functions of the 
thermodynamic driving forces fXi = ^ , the mass balance 
equations for the three components are 
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with mobility coefficients given by 
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The form of Eq. {Tj) allows different values of the bare 
trace diffusion coefficients Di{(j>) for the different compo- 
nents i and satisfies the constraint ci + C2 + C3 = 1. 

Numerical aspects of phase-field modeling have been 
discussed in |3 including the conditions for spatial grid 
resolution in comparison to the interface thickness. The 
evolution Eqs. l(2Jl and © are solved using a finite ele- 
ment method with a semi-implicit time discretization on 
a nonuniform adaptive mesh having the highest order of 
spatial resolution in the vicinity of the solid-liquid inter- 
face where the gradients of the phase field and of the con- 
centrations reach maximal values. The adaptive grid re- 
finement criterion has been defined ensuring a minimum 
of 7-10 grid points in regions of the diffuse phase bound- 
ary. To verify the origin of interface instability, simu- 
lations with different time and spatial resolution have 
been performed leading to the same growth morphology. 
Hence, it has been shown that the interfacial instabilities 
are of MuUins-Sekerka-type and not numerically induced. 

The solidus and liquidus lines of the binary Ni-Cu and 
of the Ni-Cr phase diagram can be constructed from 
the free energy, Eq. using the melting tempera- 
tures TNi = 1728 K, Tcn = 1358 K, Tcr = 1465 K, 
the latent heats Lni = 2350 J/cm^, Lcu = 1728 J/cm^, 
Lcr = 1493 J/cm"^ and the molar volume Vm — 7.42 cm'^. 
The values of the melting temperature and of the latent 
heat for Cr are adjustable parameters in order to recover 
the actual binary phase diagram in the given region of 
concentrations. The above values lead to a partition coef- 
ficient /cg — 0.843, to a liquidus slope — —3.27 K/at.% 
and to the undercooling ATq — 24.4 K for the binary 
Ni6oCu4o system. Similarly, we obtain ke = 0.905, 
TOe -2.08 K/at.% and ATq = 8.7 K for the binary 
Ni6oCr4o system on the corresponding equilibrium phase 
diagram. 

To accentuate the effect of solidification interval, we 
assume that surface properties (surface energy density cr 
and kinetic coefficient ^) do not depend on alloy com- 
position and have the values a = jqT = 0.37 J/m'^ 
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and /i = 3.3 mm/(s-K) The anisotropy of the in- 
terface properties plays an important role in the selec- 
tion of the operating state during dendritic growth. In 
this study, we use the values obtained from molecular- 
dynamics simulations for the anisotropics in pure Ni, 
Eqs. I0J and ©, = 0.023 for surface free energy density 
and Efc = 0.169 for the kinetic coefficient 0- As reported 
in , the diffusion coefhcicnts of Ni and Cu in melt are 
Dm = 3.82 X 10-^ mVs, L»cu = 3.32 x 10"^ mVs and we 
assume Dq^ — 1.5 x 10~^ m^/s. The diffusion coefficient 
in the solid phase is set equal to lO"^'^ xr? /s for all com- 
ponents. The value of the small length scale parameter 
in the entropy functional, Eqs. |^ and (O, is chosen to 
be £ = 0.1 /im. 

We conducted numerical computations for differ- 
ent alloy compositions varying from Ni6oCu36Cr4 to 
Ni6oCu4Cr36 and for fixed undercoohng conditions. We 
kept the concentration of Ni at 60at.% and adjusted the 
initial undercooling at 20 K measured from the equilib- 
rium liquidus line for a given composition. The corre- 
sponding shapes of growing crystals are shown in Fig. ^ 
At small contents of Cr, the solid phase forms a den- 
dritic structure that has a pronounced preferable growth 
direction determined by the anisotropy of surface energy 
and interface kinetics (first three images of Fig. ^| . An 
increase of Cr concentration in the melt leads to an in- 
crease of dendritic tip velocity V . This is accompanied by 
a thickening of the primary dendritic trunk. Further rise 
of Cr concentration (middle three images of Fig. ^ causes 
a morphological transition from a dendritic to a globu- 
lar growth shape at a composition Ni6oCu2oCr2o. The 
globular form of the crystals is stable after this morpho- 
logical transition for further change of alloy composition 
to Ni6oCu4Cr36. Within the globular regime of solidifica- 
tion, the anisotropy of the interface has a smaller effect 
on the crystal shape. The preferred growth directions are 
less pronounced in comparison with the dendritic shape 
for amounts of Cr less than 20at.%. To investigate the 
dependence of the spatial redistribution of the alloy com- 
ponents in the solid phase on the melt composition, we 
consider isolines corresponding to the average concentra- 
tion of nickel atoms in the solid phase. These isolines 
separate Ni depleted and Ni enriched domains (see Fig.^ 
solid lines). The character of the isolines for copper and 
chrome is qualitatively the same only having small quan- 
titative deviations. Two independent parts can be found 
in the geometry of the isolines. The first part is aligned 
in a direction parallel to the dendritic trunk remaining 
in this position after the transition to the globular mor- 
phology. The second part develops in diagonal direction 
and it disappears at the morphological transition. 

The morphological transition can quantitatively be 
measured using the deviation (5 = (5 — 5o)/S'o of area S 
covered by the solid phase from the triangular area So be- 
tween the dendritic tips and the center of the initial seed. 
Fig. 121 With the change of melt composition, the devia- 



tion 5 changes the sign from a negative value for dendritic 
morphologies to a positive value for globular shapes. The 
transitional value (5 = corresponds to the composition 
Ni6oCu2oCr2o and is hence in good agreement with the 
transition point given by the vanishing diagonal part of 
the isoline. The transition takes place in a smooth man- 
ner. On the one hand, this follows from the linear increas- 
ing of the growth velocity V of the dendritic or globular 
tip with increasing Cr concentration (see the values for 
V in Fig.^l. On the other hand, the smooth transition is 
confirmed by the smooth dependence of the deviation 5 
on the melt composition, Fig. |2 The morphology of the 
growing structure is determined by the concurrence of 
the stabilizing force due to the surface energy and of the 
destabilizing force due to the concentration gradients on 
the interface P, Q • The morphological transition can be 
explained considering the balance between the stabilizing 
and destabilizing forces. During the computed transition 
from a dendritic to globular shape for varying composi- 
tion of the ternary alloy Ni6oCu4o-a;Cr2, with < x < 40, 
the liquidus slope and the deviation of the equilibrium 
partition coefficient from unity decrease. This leads to a 
reduction of the destabilizing forces on the interface. The 
globular growth shapes are becoming stable against per- 
turbations. Although the phase-field simulations of the 
ternary system were conducted for fixed undercooling, 
we remark that a similar morphological transition from 
dendritic to globular growth shapes has been observed in 
numerical simulations of dendritic growth in binary al- 
loys ^3|> where the initial undercooling acts as a control 
parameter for the transition. In this case a morphological 
transition occurs through a coalescence of side branches 
with increasing of the undercooling. 

In summary, our numerical simulations show the oc- 
currence of a morphology transition from dendritic to 
globular growth structures in the ternary Ni-Cu-Cr sys- 
tem for a varying alloy composition. The results demon- 
strate the general potential of the phase-field approach 
in modeling multicomponent solidification in a thermo- 
dynamically consistent way. Together with previous de- 
velopments ,.5J, the work in this Letter shows that the 
diffuse interface formulation provides a uniform descrip- 
tion for modeling complex microstructure evolution in 
pure, binary and multicomponent alloy systems. 
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FIG. 1: Transition from dendritic to globular morphology for different initial melt compositions at a constant initial undercooling 

of 20 K. Shaded regions correspond to the solid phase, i.e. (j) ^ 0.5, and solid lines represent the isolines of average concentration 
of Ni in solid. The values for the tip velocities V and for the area estimations 6 are stated in each image. 
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FIG. 2: Simulated area S covered by the crystal in relation 
to a normalized triangular area So for different melt compo- 
sitions. The 5 = line marks the transition from dendritic to 
globular growth morphologies. 
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